Time reversal, fermion doubling, and the masses of lattice Dirac fermions in three 

dimensions 
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Motivated by recent examples of three-dimensional lattice Hamiltonians with massless Dirac 
fermions in their (bulk) spectrum, I revisit the problem of fermion doubling on bipartite lattices. The 
number of components of the Dirac fermion in a time-reversal and parity invariant d-dimensional 
lattice system is determined by the minimal representation of the Clifford algebra of d + 1 Her- 
mitian Dirac matrices that allows a construction of the time-reversal operator with the square of 
unity, and it equals 2 d for d = 2, 3. Possible mass-terms for (spinless) Dirac fermions are listed and 
discussed. In three dimensions there are altogether eight independent masses, out of which four 
are even, and four are odd under time reversal. A specific violation of time-reversal symmetry that 
leads to (minimal) four-component massless Dirac fermion in three dimensions at low energies is 
constructed. 



I. INTRODUCTION 

It is well-known that the tight-binding Hamiltonian 
for spinless fermions on a honeycomb lattice exhibits two 
inequivalent zero-energy points in the Brillouin zone, in 
vicinity of which the excitation energy has an approx- 
imate linear momentum dependence Honeycomb 
lattice being bipartite, this implies that the low-energy 
excitations behave like four-component Dirac fermions. 
Another well-known two-dimensional example of four- 
component Dirac fermions arising on a square lattice is 
the 7r-fl.ux Hamiltonian Q , in which there is a single Dirac 
point, but the unit cell consists of four lattice sites. Of 
course, this is precisely how it has to be according to the 
Nielsen- Ninomiya (NN) theorem 3] , which forbids the ap- 
pearance of two-component (Weyl) fermions on a lattice 
under a certain set of conditions, which, importantly, in- 
cludes the time-reversal invariance of the Hamiltonian. 

Several ways of circumventing the NN "no-go" theo- 
rem do exist, however. It was pointed out Q that it may 
be possible to find a two-component massless fermion 
on a two-dimensional lattice at the expense of having 
additional low-energy excitations in the system, such as 
flat bands, for example. Furthermore, breaking the time- 
reversal and then fine-tuning the hopping on the honey- 
comb lattice can also lead to a single Weyl fermion at 
low-energies ji 5] • Since neutrinos are such Weyl fermions, 
a possibility of their lattice realization has over the years 
attracted a considerable attention in the lattice field the- 
ory community 6]. In condensed matter physics Weyl 
fermions would represent a new class of quasiparticle 
excitations with interesting properties, presently under 
scrutiny in the context of topological band insulators, 
where they appear generically as surface modes. Q 

Recently, interesting generalizations of both 7r-flux and 
graphene's lattice Hamiltonians to three dimensions (3D) 
have appeared [8| , both displaying (for spinless fermions) 
eight-component Dirac fermions as their low-energy ex- 
citations. Since the number of order parameters, or the 
mass terms, in 3D is larger than in 2D, it seems there 
could be more possibilities for a reduction of the dimen- 



sion of the low-energy Dirac fermion by fine tuning the 
masses a la Haldane. This issue is investigated here, to 
find that generally the "smallest" Dirac fermion that can 
be found in the spectrum by breaking the time-reversal in 
3D still has four components. Along the way, it is shown 
that the minimal dimension of 2 d of the Dirac fermion in 
a d-dimensional time-reversal-symmetric lattice Hamil- 
tonian can also be understood as a purely algebraic con- 
straint: a smaller representation of Dirac matrices does 
not allow a construction of a time-reversal operator for 
the Dirac fermion that would square to +1, as it must 
for the spinless fermions hopping on a lattice. 

A 2 d -component low-energy massless Dirac Hamilto- 
nian in d-dimensions exhibits 0(d + 1) "chiral" symme- 
try that may be used to classify possible (order param- 
eter) mass terms. In a representation-independent way, 
I show that in general there are d + 1 time-reversal in- 
variant mass terms, which form a vector under the chiral 
rotations. In addition, there are also one (in d = 2) and 
four (in d = 3) masses that break time-reversal invari- 
ance, and which are scalar and vector under the chiral 
symmetry, respectively. A possible physical meaning of 
these masses on the 2D honeycomb [l|, [H, l9j-[T]| and the 
2D [12J, |l3[ and 3D [8[ 7r-flux lattices has been discussed 
at length elsewhere. The commutation relations between 
the mass terms in 3D allow them to be fine-tuned so as 
to reduce the number of components of the low-energy 
Dirac excitation, albeit only to four. 



II. FERMION DOUBLING: ALGEBRAIC 
DERIVATION 



Assume that all low-energy excitations of a lattice 
Hamiltonian for spinless electrons hopping in d spatial 
dimensions may be described by an effective Dirac Hamil- 
tonian 



H = J2 



CtiPi 



(1) 
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where pi are the usual momentum operators, a, are 
Hermitian Dirac matrices satisfying the Clifford algebra 
{cti,ctj} = 25ij, of yet to be determined dimensionality. 
The Fermi velocity in the above expression is set to unity. 

Let us deduce the dimension of a-matrices in the effec- 
tive Dirac Hamiltonian in 3D [d — 3), under two assump- 
tions. First, assume that inversion symmetry is obeyed 
by the underlying lattice Hamiltonian. Since it then also 
needs to be respected within the low-energy Dirac sector, 
it means that there exists a fourth, parity matrix, (3, such 
that 

{/W = 0- (2) 

Note that this is automatically the case for the bipar- 
tite lattice Hamiltonians. Second, let us assume the 
time-reversal symmetry. Since the lattice Hamiltonian 
describes the real-space motion of spinless particles, in 
the basis of lattice sites the antilinear operator It for 
the time-reversal is simply that of complex conjugation, 
which then satisfies if — 1. Since, on the other hand, 
the value of If is basis independent [14| , it is the same in 
the basis of Bloch states. This means that for a time- 
reversal lattice Hamiltonian there must exist an operator 
I t in the low-energy Dirac sector which squares to unity. 
I show next that this immediately implies that the Dirac 
Hamiltonian must be at least eight-dimensional, i. e. the 
usual "fermion doubling" . 

The mathematical basis behind the above assertion is 
the Okubo's theory of real representation of Clifford alge- 
bras [lH . Here I reproduce only the part relevant for the 
present purposes. Define a Clifford algebra C(p, q) otp+q 
mutually anticommuting objects, p (q) of which square 
to +f (—1). The smallest purely real representation of 
C(p,q) is then 2 9 -dimensional, and it exists when p = q 
or p = q+1 (mod 8). Okubo's theory, of course, contains 
a great deal more than this corollary. An important im- 
plication of this result is that a 2™-dimensional complex 
representation of the Clifford algebra C(2n + 1,0), where 
all elements square to +1, as in the Dirac Hamiltonian 
in Eq. 1, can always be chosen so that n matrices are 
imaginary, and the remaining n + 1 real [16j . It is easy 
to check that this is indeed the case for n = 1 and n = 2, 
for example, by a direct construction. 

Returning to the 3D Dirac Hamiltonian, let us con- 
struct the time-reversal operator for the low-energy exci- 
tations. The matrices {oti, Oa, 03, /3} form the Clifford 
algebra C(4, 0), and the minimal complex representa- 
tion, as well-known, is four-dimensional. We can form 
maximally five mutually anticommuting matrices by con- 
structing the fifth (Hermitian) matrix aia2a3/3, which 
together with the first four then forms the Clifford al- 
gebra C(5,0). Okubo's general theory, or a simple in- 
spection in this case, implies then that out of these five 
four-dimensional matrices we can have at most two as 
purely imaginary. So let us chose a.\ and /3 real, and a.i 
and «3 imaginary, for example. All other choices are, 
of course, related to this one by a unitary transforma- 
tion. The time-reversal operator is It = U K, where U is 



unitary and K stands for complex conjugation. It then 
follows that for our choice the unitary part must satisfy 

{U,ai} = [U,a 2 ] = [U,a 3 ]=0, (3) 

which has two solutions: U = or U = e 1 ^ 0.20,3 P, 

with <f> arbitrary. Demanding that a spatial symmetry 
such as parity commutes with the time-reversal in the 
low-energy Dirac sector, as it does in the microscopic 
lattice Hamiltonian, selects the latter solution. Irrespec- 
tively of this condition it follows that 

It = -1, (4) 

which stands in contradiction to the lattice origin of 
the Dirac Hamiltonian. One concludes that a four- 
dimensional Dirac fermion cannot represent all the low- 
energy excitation of a time and space-reversal invariant 
Hamiltonian on a 3D lattice. 

A (reducible) eight-component representation however, 
is compatible with the time-reversal symmetry: we can 
choose all three a-matrices to be imaginary, for example, 
so that the unitary operator U commutes with all three. 
This implies that U — e^, or U — a 102 03. In ei- 
ther case it is now if — +1. Commutation of the parity 
matrix /3 and the time-reversal in the low-energy Dirac 
sector implies further that U = e 1 f Choosing the phase 
= then leaves a simple It = K. 

Similar reasoning can straightforwardly be applied to 
d = 2 as well. Assuming a±, 02 and (3 to be two- 
dimensional matrices, implies that one of them can be 
chosen imaginary, and the other two real. Choosing /3 
imaginary makes it obvious that the time-reversal oper- 
ator that would also commute with parity does not even 
exist. Going to a four-dimensional representation, on the 
other hand, allows one to chose a\ and 02 imaginary and 
/3 real, so that time-reversal operator becomes complex 
conjugation, the same as in 3D. 

III. MASSES AND SYMMETRIES 

Let us now consider possible mass-terms that upon ad- 
dition to the Dirac Hamiltonian in Eq. 1 would gap its 
spectrum. On a lattice such masses would break some 
of the symmetry and represent certain order parameter. 
Before considering three dimensions, let me review the 
familiar example of 2D first. 

A. Two dimensions 

If d = 2, we can chose the two four-dimensional ma- 
trices a\ and a 2 to be imaginary, as discussed above. 
This leaves three other linearly independent real Her- 
mitian four-dimensional matrices {/3i, $2, which to- 
gether with a-matrices form the maximal Clifford algebra 
C(5,0) of dimension four. One of the /3-matrices may 
be chosen to represent parity operator, and of course, 
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ft = aic^ftft- The time-reversal operator is then just 
the complex conjugation, and the most general time- 
reversal invariant mass-term is to ■ Vr, where 

^ = (ft, ft, ft) (5) 

is a vector under the 0(3) chiral symmetry of the Dirac 
Hamiltonian, generated by iftft, for i =fi j. The masses 
to = (mi, W2, TI3) are real. 

Multiplying an odd number of ftmatrices produces an- 
other matrix that anticommutes with a\ and 012, and 
thus leads to a new linearly independent mass term. For 
d = 2 there is only one such product, 

Si = iftftft, (6) 

where the imaginary unit is necessary to make the matrix 
Hermitian. This mass is evidently a scalar under 0(3) 
symmetry, but odd under time-reversal. On the honey- 
comb lattice Vr represents a charge density wave, and 
two components of a Kekule bond density wave [llj. Si, 
on the other hand, represents the anomalous quantum 
Hall stated]. 

B. Three dimensions 

In d = 3, we chose again the three eight-dimensional 
Hermitian matrices a*, i — 1, 2, 3 to be imaginary, to have 
simply I t — K. Now, however, there are two, related but 
distinct ways to complete this set to the maximal number 
of seven mutually anticommuting matrices. First, we can 
chose four real matrices ft, i — 1,2, 3, 4, which together 
with a-matrices form O (7, 0) . Any linear combination 
to • Vr of these would yield a time- reversal invariant mass 
term, with 

V R = 081, ft, ft, ft), (?) 

again as a real vector under 0(4) symmetry of the Dirac 
Hamiltonian. This is analogous to the situation in d = 2. 
The novelty in d — 3 is that there are now four dis- 
tinct three-products of /3-matrices, which are also possi- 
ble mass terms. Consider a 0(4) vector 

^ = (ft,ft,ft,ft), (8) 

with ft = iftftftftft. Vi is clearly odd under time- 
reversal, due to the imaginary unit necessary to render 
the matrices Hermitian. 

A concrete realization of the above algebra on a cu- 
bic lattice with a 7r-flux through a plaquette was recently 
proposed and discussed by Hosur, Ryu, and Vishwanath 
|8|]. Time-reversal-even matrices in their example rep- 
resent the charge-density and bond-density wave order 
parameters, and the time-reversal-odd order parameters 
are three anomalous Hall states and a "chiral" topologi- 
cal insulator. 



IV. FOUR-COMPONENT MASSLESS FERMION 
IN 3D: A CONSTRUCTION 

Finally, let us now investigate the possibility of having 
a Dirac massless fermion with fewer than eight compo- 
nents in the spectrum by a judicious breaking of time- 
reversal symmetry. I consider only d — 3 in detail, al- 
though the result for d = 2, which has been pointed out 
by Haldane long ago, would be analogous. 

In d — 3 the /3-matrices obey the commutation rela- 
tions: 

{ft,ft} = {ft,ft} = 2%, (9) 

{ft,ft} = [ft,ft ¥ ;]=0. (10) 

Consider then the most general massive Dirac Hamilto- 
nian in d=3: 

3 

H = ^a lPl + rh-V R + m-Vi. (11) 

i=l 

Without loss of generality we can utilize the 0(4) chi- 
ral invariance and chose Vr = (ft, 0,0,0), so that to = 
(to, 0, 0, 0). Squaring the Dirac Hamiltonian we then find 

4 

H 2 =p 2 +TO 2 +TO 2 +TO^TO l {ft,ft} (12) 

1=2 

where to 2 = X^i=i Consider the three matrices 

7* = ^{ft,ft} = ftft, (13) 
i = 2,3,4, that appear in the last expression. Evidently, 
{7^} = 2% (14) 

and they themselves close the Clifford algebra O(3,0). 
Since the two inequivalent irreducible representations of 
0(3, 0) are two dimensional, the three 7-matrices can be 
brought into a block-diagonal form, with the blocks as 
Pauli matrices. This implies that there exists a unitary 
transformation that puts the (square of) Hamiltonian 
into the form 

H 2 = p 2 +m 2 +m 2 +2m(m 2 2 +ml+rhl) 1/2 (I 4 ®<j 3 ), (15) 

where £4 is a four-dimensional unit matrix. The energy 
spectrum is therefore 

E 2 ± =p 2 + {m± (ml +m\ + m\) 1/2 ) 2 + to 2 , (16) 

with each sign being four times degenerate. 

Since choosing the mass Vr breaks the 0(4) symme- 
try, the obtained spectrum depends qualitatively on the 
"direction" of the vector Vj. Choosing Vi parallel to Vr 
adds the corresponding masses in squares and only in- 
creases the mass gap. Taking Vj to be orthogonal to Vr, 
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on the other hand, increases the mass gap for one half, 
and decreases it for the other half of the components of 
the original eight-component Dirac fermion. In particu- 
lar, for rhi = and 

m 2 = fh\ + m| + m| (17) 

the mass-gap for four of the components collapses to zero, 
and one finds a four-component massless Dirac fermion 
as a sole low-energy excitation. It is obvious from the 
above derivation, however, that no combination of eight 
mass-terms could leave a two-component Weyl fermion 
at low energies in d = 3. This is possible only in d = 2, 
where the Weyl fermion results from halving the original 
four-component Dirac fermionQ, in precise analogy to 
the above derivation in 3D. 



V. CONCLUSION AND DISCUSSION 

I studied massless Dirac fermions as the low-energy 
excitations of lattice Hamiltonians in two and three di- 
mensions in a general, representation-independent way. 
The existence of the time-reversal operator with a trivial 
square in the low-energy sector is shown to imply that the 
Dirac fermion must have at least four and eight compo- 
nents, in 2D and 3D, respectively. In 3D there are eight 
possible mass terms for (spinless) Dirac fermions, four 
of which violate time-reversal symmetry. A combination 
of time-reversal-even and time-reversal-odd masses that 
leaves a four-component massless Dirac fermion in the 
spectrum in 3D was constructed. 

Recently, an interesting variation of the 7r-flux Hamil- 
tonian in 2D was proposed [l3T | which seemingly has two 
Weyl fermions with different velocities as its low-energy 
spectrum, although it is perfectly time-reversal and par- 
ity invariant. It is interesting to see why this does not 
contradict our assertion that the minimal dimension of a 
Dirac fermion in d dimensions is 2 d . The Hamiltonian of 
the ref. 13 can be written as 

H = (ai + 5ihh)pi + (a 2 + Sipxps)p2 (18) 

in the notation of sec. Ill a, with 5 as a real parameter. 
The energy spectrum is then ±(1 ± S)p, and choosing on 
imaginary and /3j real, the time-reversal operator which 



commutes with the Hamiltonian is I t = K, in accord with 
the lattice origin of the Hamiltonian. The contradiction 
with our assertion is avoided by the fact that H cannot 
be transformed into a block-diagonal form that would 
describe an orthogonal sum of two Weyl Hamiltonians. 
This can be proven by reductio ad absurdum. If it could 
be so block-diagonalized, it would have to be equivalent 
to 

H d = (l + 5)(piai +p 2 o- 2 ) © (1 - S)( Pl a! + p 2 o- 2 ), (19) 

and therefore there would exist a traceless matrix that 
squares to unity, and that commutes with H for any mo- 
mentum, which would be equivalent to 03 ® oq. But, 
such a matrix obviously does not exist: the only matrix 
that commutes with a±, a 2 , and ifi\fi 2 is «/?i/?2 itself, 
which however fails to commute with i/?i/?3. The exis- 
tence of matrix f3\ which anticommutes with the Hamil- 
tonian in Eq. 18, on the other hand, implies that it can 
be transformed into a block off-diagonal form. This is a 
nice example of how the presence of additional degrees 
of freedom at arbitrary low energy only seems to invali- 
date the NN theorem. Note that putting the parameter 
5 = 1 in the above Hamiltonian makes two of the bands 
flat and with zero energy, somewhat similar to the early 
example of Dagotto, Fradkin, and Moreo [4|. 

Several extensions of the present work can be envis- 
aged. One can include the particle's spin and/or consider 
off-diagonal pairing terms in 3D to determine all possible 
Dirac mass terms, similarly to what has been done in 2D 
[l7j | . Extensions to higher dimensions, albeit somewhat 
academic from the point of view of a condensed matter 
physicist, may be interesting from a particle physics or a 
methodological perspective. Finally, it could be interest- 
ing to investigate the structure of topological defects in 
the mass-terms in 3D, as several analogies with the more 
familiar defects in 2D [l8[ are known to exist. [ijj H3] 
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